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Abstract. Heavy- ion reactions provide a unique means to investigate the equation 
of state (EOS) of neutron-rich nuclear matter, especially the density dependence of the 
nuclear symmetry energy Esym if) ■ The latter plays an important role in understanding 
many key issues in both nuclear physics and astrophysics. Recent analyses of heavy-ion 
reactions have already put a stringent constraint on the Esym (p) around the saturation 
density. This subsequently allowed us to constrain significantly the radii and cooling 
mechanisms of neutron stars as well as the possible changing rate of the gravitational 
constant G. 



PACS numbers: 26.60.-|-c,24.10.-i 
1. Introduction 

Within the parabolic approximation which has been verified by all many-body theories 
to date, the EOS of isospin asymmetric nuclear matter can be written as E{p, 5) = 
E{p,6 = 0) + i?sym(p)5^ + where 6 = {pn - Pp)/{pp + Pn) is the isospin 

asymmetry and Esym{p) is the density-dependent nuclear symmetry energy. The latter 
is very important for understanding many interesting astrophysical problems [H [2] , the 
structure of radioactive nuclei Oil], and the dynamics of heavy- ion reactions O [6l [TJ [8] . 



2. Constraining the density dependence of the nuclear symmetry energy in 
terrestrial nuclear laboratories 

Considerable progress has been made recently in determining the density dependence of 
the nuclear symmetry energy at sub-normal densities in terrestrial nuclear laboratories. 
Shown in Fig. [1] are several typical theoretical model predictions P [THl [H] compared 
with the constraints obtained from model analyses of experimental data (those labeled 



Constraining properties of neutron stars with heavy-ion reactions in terrestrial laboratories2 




Figure 1. Left panel:Density dependence of the nuclear symmetry energy using 
the MDI interaction with x = and x = —1 and other many-body theories 
predictions (taken from |22j). Right paneh symmetry energies obtained from 
21 sets of Skyrme interactions and the MDI interaction with x = —1 and 
X = [19]. 



a; = 0, a; = — 1 and FSU-Gold). The constraints labeled x = and a; = — 1 were 
extracted recently from studying isospin diffusion in the reaction of ^^''Sn +^^^Sn at 
Ebeam/A = 50 MeV within a transport model [121 [131 [HI [15] using the MDI interaction 
|16j . For this particular reaction the maximum density reached is about 1.2po- The 
constraints are thus only valid below this density. Moreover, it was shown that the 
neutron-skin thickness in 208p]^ calculated within the Hartree-Fock approach using 
the same underlying Skyrme interactions as the ones labeled x = and x = — 1 is 
consistent with the available experimental data [T71 [181 [IH]- The symmetry energy 
labeled as FSU-Gold was calculated within a Relativistic Mean Field Model (RMF) 
using a parameter set such that it reproduces both the giant monopole resonance in 
^°Zr and ^^'^Pb, and the isovector giant dipole resonance of ^ospi-, j2Qj_ These results 
all together represent the best phenomenological constraints available on the symmetry 
energy at sub-normal densities. The various predictions at supra-normal densities still 
diverge widely. Hopefully, nuclear reactions with high energy radioactive beams will 
allow us to pin down the high density behavior soon using several probes predicted 
recently [2T] . 

The available constraints on the symmetry energy limit the nuclear effective 
interactions in nuclear matter. This can be seen by comparing them with the symmetry 
energies obtained from Skyrme effective interactions [I9]. The right panel of Fig. [1] 
displays the density dependence of Esy^{p) for 21 sets of Skyrme interaction parameters 
that are currently used widely in nuclear structure studies. Surprisingly, most of these 
effective interactions lead to symmetry energies rather inconsistent with the constraints 
discussed above. 



Constraining properties of neutron stars with heavy-ion reactions in terrestrial laboratories?) 



3. Constraining the radii and cooling mechanisms of neutron stars 



While the maximum mass of neutron stars is mainly determined by the incompressibility 
of symmetric nuclear matter, their radii are primarily determined by the isospin 
asymmetric pressure that is proportional to the slope of the symmetry energy E'^^^{p). 
For the simplest case of a neutron-proton-electron {npe) matter in neutron stars at (3 
equilibrium, the pressure is given by P(p, 6) = Po(p) + -Pasy(p, = (|f )^ + iPePe 
= E'{p,6 = 0) + E'^y^{p)6'^ + — 6)pEsy^{p), where the electron density is 
Pe = |(1 — S)p and the chemical potential is pe = Pn — Pp = 45-Esym(p)- The equilibrium 
value of S is determined by the chemical equilibrium and charge neutrality conditions, 
i.e., 6 = 1 — 2xp with Xp ~ 0.048 [-Esym(p)/-Ssym(Po)]'^ (p/Po)(l ~ 2xp)^. Because of the 
large 6 value in neutron stars, the electron degenerate pressure is small. Moreover, 
the isospin symmetric contribution to the pressure is also very small around normal 
nuclear matter density as E'{po,6 = 0) = 0. The pressure is thus dominated by the 
term proportional to the slope of the symmetry energy. Since neutron star radii are 
determined by the pressure at moderate densities where the proton content of matter is 
small, they are very sensitive to the slope of the symmetry energy near and just above 
Pq. In particular, a stiffer symmetry energy is expected to lead to a larger neutron star 
radius. 
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Figure 2. Left panel: mass and proton fraction of neutron star and the 
symmetry energy as functions of density. Right panel: correlations in neutron 
star masses and radii [18]. All results are for spherically-symmetric, non- 
rotating, non-magnetized neutron stars consisting of npep matter at zero 
temperature. 



The dependence of some basic properties of neutron stars on the nuclear EOS is 
shown in the left panel of Fig. [3l The top and middle panels give, respectively, the 
mass of neutron star and its proton fraction Xp, calculated from the MDl interaction 
with a; = 0, —1, and —2, as functions of its central density. The symmetry energies 
obtained from these interactions are shown in the lower panel together with that from 
the AV18-|-5u-|-UIX* interaction of Akmal et al. (APR) [23]. It is interesting to see 
that up to about 5po the symmetry energy predicted by APR agrees very well with that 
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from the MDI interaction with x = 0. For Xp below 0.14 [Ij, the direct URCA process 
for fast coohng of proto-neutron stars does not proceed because energy and momentum 
conservation cannot be simultaneously satisfied. For EOSs from the MDI interaction 
with X = —1 and x = —2, the condition for direct URCA process is fulfilled for nearly all 
neutron stars above 1 Mq. For the EOS from x = 0, the minimum density for the direct 
URCA process is indicated by the vertical dotted line, and the corresponding minimum 
neutron star mass is indicated by the horizontal dotted line. While this indicates that 
neutron stars with masses above 1.39 Mq have a central density above the threshold for 
the direct Urea process, this conclusion may be modified by the presence of terms in 
the symmetry energy which are quartic in the isospin asymmetry at high density [2l]. 

The relations between neutron star masses and radii for above EOSs are given 
in the right panel of Fig. [3l Also given are the constraints due to causality as well 
as the mass-radius relations from estimates of the crustal fraction of the moment of 
inertia (A/// = 0.014) in the Vela pulsar [25] and from the redshift measurement from 
Ref. [2^. Allowed equations of state should lie to the right of the causality line and 
also cross the other two lines. The hatched regions are inferred limits on the radius 
and the radiation radius (the value of the radius which observed by an observer at 
infinity) defined as Roo = R/ \/l ~ 2GM/ Rc^ for a 1.4 Mq neutron star. It is seen 
that the symmetry energy affects strongly the radius of a neutron star but only slightly 
its maximum mass [H [27]. These analyses have led to the conclusion that only radii 
between 11.5 and 13.6 km (or radiation radii between 14.4 and 16.3 km) are consistent 
with the EOSs from the MDI interaction with x = and x = —1 and thus with the 
terrestrial nuclear laboratory data. The observational determination of the neutron star 
radius from the measured spectral fluxes relies on a numerical model of the neutron star 
atmosphere and uses as inputs the composition of the atmosphere, a measurement of the 
distance, the column density of x-ray absorbing material, and the surface gravitational 
redshift. Since many of these quantities are difficult to measure, only a paucity of radius 
measurements are available. Nevertheless, it is interesting to note that the calculations 
shown in Fig[3] are consistent with recent observations p8] including the very rapidly 
rotating pulsar discovered recently [29j. 



4. Constraining the changing rate of the gravitational constant G 

Testing the constancy of the gravitational constant G has been a longstanding 
fundamental question in natural science. As first suggested by Jofre, Reisenegger and 
Fernandez [30], Dirac's hypothesis |3T]| of a decreasing gravitational constant G with 
time due to the expansion of the Universe would induce changes in the composition 
of neutron stars, causing dissipation and internal heating. Eventually, neutron stars 
reach their quasi-stationary states where cooling due to neutrino and photon emissions 
balances the internal heating. As shown in ref. [SDj the stationary surface temperature 
is directly related to the relative changing rate of G via = V\^\ , where the 
function X> is a quantity depending only on the stellar model and the equation of state. 
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Figure 3. Left: Neutron star stationary surface temperature for stellar models 
satisfying the mass constraint by van Straten et al. [33] . The solid lines are the 
predictions versus the stellar radius for the considered neutron star sequences. 
Dashed lines correspond to the 68% and 90% confidence contours of the black- 
body fit of Kargaltsev et al. [34J. The value of \G/G\ = 4 x IQ-^^y^-i 
chosen so that predictions from the x = EOS are just above the observational 
constraints. Right: Same as left but assuming \G/G\ = 8 x 10~^^yr~^ . 



The correlation of surface temperatures and radii of some old neutron stars may thus 
carry useful information about the changing rate of G. Using the constrained symmetry 
energy with x = and x = —1 shown in Fig{Tl within the gravitochemical heating 
formalism, as shown in Fig|3]we obtained an upper limit of the relative changing rate 
of \G/G\ < 4 X 10~^^yr~^. This is the best available estimate in the literature [32] ■ 
For a comparison, results with the EOS from the recent Dirac-Brueckner-Hartree-Fock 
(DBHF) [35] calculations using the Bonn B potential are also shown. The Bonn B 
potential gives roughly the same result on the stationary surface temperature but slightly 
larger radius compared to the x = case. 

5. Summary 

In summary, the symmetry energy of neutron-rich matter is fundamentally important 
for both nuclear physics and astrophysics. Available data from heavy-ion reactions in 
terrestrial laboratories allowed us to constrain the symmetry energy at sub-saturation 
densities. This has led to stringent constraints on not only the nuclear effective 
interactions but also the cooling mechanisms, radii of neutron stars and the changing 
rate of the gravitational constant G. 
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